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TeV and Superheavy Mass-Scale Particles from Supersymmetric Topological Defects, 
the Extragalactic 7-ray Background, and the Highest Energy Cosmic Rays 
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Cosmic topological defects in a wide class of supersymmetric theories can simultaneously be sources 
of higgs particles of mass comparable to the supersymmetry breaking scale ~ TeV, as well as su- 
perheavy gauge bosons of mass ~ r\, where n (3> 1 TeV) denotes the associated gauge symmetry 
breaking scale. For cosmic strings with r\ ~ 10 14 GeV, the higgs decay can account for the extra- 
galactic diffuse 7-ray background above ~ 10 GeV, while the gauge boson decay may explain the 
highest energy cosmic ray flux above ~ 10 11 GeV, provided that particle production is the dominant 
energy loss mechanism for cosmic strings, as recent simulations suggest. By the same token, cosmic 
strings with 77 much above ~ 10 14 GeV will be ruled out. 

PACS numbers: 98.80.Cq, 12.60.Jv, 98.70.Sa, 98.70.Vc 



In a wide class of supersymmetric (SUSY) unified gauge 
theories, including some versions of effective theories de- 
rived from superstrings, certain phase transitions can oc- 
cur at a temperature comparable to the "soft" super- 
symmetry breaking scale < 1 TeV, even though the asso- 
ciated gauge symmetry breaking scale itself may be much 
larger jl]]^] . This occurs rather generically in SUSY theo- 
ries as a consequence of existence [[| of 'directions' along 
which the effective potential V of the relevant scalar field 
is almost flat, i.e., the curvature \y"\^l 2 of the poten- 
tial is much smaller than the vacuum expectation value 
(VEV) 77 of the scalar field out to field values ~ 77. An 
"almost flat potential" for a (complex) higgs field $ after 
supersymmetry breaking generally has the form |l],|],(||5) 
V = V - m 2 2 + X n mpf n ^ 2n +\ where = |$| 

and mpi = (87TG) 1 / 2 ps 2.4 x 10 18 GeV. The (j) 2 term 
arises from 'soft' supersymmetry breaking, so the mass 
scale m s is typically < 1 TeV. The higher order (non- 
renormalizable) terms would arise from 'integrating out' 
particles of Planck mass scales in a 'higher' theory such 
as superstring theory. The flatness of the potential is due 
to absence of the Xcj) 4 term familiar in non-SUSY theo- 
ries. Depending on the strengths of the couplings A„, the 
minimum of the potential, i.e., the VEV 77, can lie any- 
where in the range - 10 9 GeV to - M GU T ~ 10 16 GeV, 
the grand unification (GUT) scale The 'height' of the 
potential is Vq ~ rj 2 m 2 . 



For temperatures T in the range m s <C T < U c 
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finite temperature corrections to the potential can hold 
the <& field at 6 = until T falls below m s , at which 



a phase transition occurs taking <f> to 77. Note that al- 
though the higgs scalars in these theories are 'light' with 
mass m,p ~ \V"\ 1 / 2 ~ m s < 1 TeV, the associated gauge 
bosons have the 'usual' mass ~ 77, which can, in particu- 
lar, be ~ Mqut, if the higgs under consideration breaks 
the GUT symmetry. Some cosmological consequences of 
theories with almost flat potentials (hereafter simply re- 
ferred to as "flat SUSY" theories) have been considered 
inRefs. |g§. 

In this Letter, we point out that cosmic topological 
defects (TDs) @ such as magnetic monopoles and cosmic 
strings associated with phase transitions in flat SUSY 
theories |^|| can, through their collapse, annihilation, 
or other processes, be sources of higgs bosons of mass 
~ TeV, as well as of gauge bosons of superheavy mass 
scale ~ 77 ^> 1 TeV, and that the decay products of both 
these kinds of particles may be observable in the Universe 
today. 

Production of extremely energetic photons, nucleons 
and neutrinos through decay of massive "X" particles 
(typically of GUT-scale mass ~ 10 16 GeV) originating 
from TDs §-@, is a subject of much current interest 
as a possible explanation [|l3] |l5| of the highest energy 
cosmic ray (HECR) events at energies > 10 11 GeV p6| . 
These TDs have usually been considered within the con- 
text of the standard non-SUSY (and non-flat) quartic 
GUT symmetry-breaking higgs potential 0, for which 
the relevant phase transition occurs at a temperature 
T ~ 77 ~ 10 16 GeV. In this case, the associated GUT 
gauge bosons as well as higgs bosons, and consequently 
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the X particles 'constituting' the TDs, all have masses 
of order the GUT scale VEV 77 ~ 10 16 GeV. In con- 
trast, the new feature in flat SUSY theories is that the 
X particles produced by the same TD processes |] [lj] 
would now be higgs of mass ~ TeV as well as super- 
heavy gauge bosons of mass ~ 77 ^> 1 TeV. We show 
that decay of these TD-produced TeV mass-scale higgs 
may contribute significantly to the extragalactic diffuse 
7-ray background (EDGRB) Q above a few GeV (which 
seems to be difficult to explain otherwise in terms of emis- 
sions from astrophysical objects; see below), while the 
superheavy gauge bosons could be a source of the HECR 
particles. In particular, we show that cosmic strings in 
flat SUSY theories with 77 ~ 10 14 GeV may simultane- 
ously explain both EDGRB above a few GeV and HECR, 
if X particle production (rather than gravitational radi- 
ation emission) is the dominant energy loss mechanism 
for cosmic strings — a possibility recently suggested in 
Ref. j[o|. By the same token, cosmic strings with 77 
much larger than 10 14 GeV (and hence GUT-scale cos- 
mic strings with 77 ~ 10 16 GeV) in flat SUSY theories 
overproduce both EDGRB and HECR, and are, there- 
fore, ruled out. (In this case, cosmic strings with 77 much 
larger than 10 14 GeV in non-SUSY theories are also ruled 
out because they overproduce HECR). 

Note that although the higgs in flat SUSY theories are 
'light', they are not expected to be produced in acceler- 
ators operating at energies well below the energy scale 
~ rj because their coupling to minimal supersymmetric 
standard model (MSSM) particles, for example, is ex- 
pected to be suppressed by a factor ~ m^/rj <C 1. Thus, 
TDs may indeed be the only source of these higgs in the 
present day Universe. 

The rate of X particle production per unit volume at 
a time t in the matter dominated epoch from a sys- 
tem of TDs can be generally written in the form ||] 
dnx/dt = (Qo/mx) {t/ta)~ A+p , where to denotes the 
present epoch, and Qo is the rate of total energy re- 
leased per unit volume in the form of X particles (higgs 
plus gauge bosons) in the present epoch. (We use 
natural units, K = c = 1, throughout.) On general 
grounds we expect that the total energy released will be 
roughly equipartitioned between the higgs (Qo.0) an d the 
gauge boson (Qo,v) modes, and so we will assume that 
Qo,<t> ~ Qo.v ~ (l/2)<3o- The dimensionless parameter 
p is in general different for different systems of TDs || . 
Here we consider the case p = 1, which is representative 
of a large class of TD processes including those involv- 
ing cosmic strings and magnetic monopoles ^ [HJ . The 
case in which the X particles are of heavy mass scales 
~ 0(77) has been considered earlier (||],[ll] |H|- Here 
we consider the effects of the TeV mass-scale higgs X 
particles, which we shall assume to be non-relativistic. 
These higgs would decay on a time scale r ~ 

6.6^ 1 (^/lO 16 GeV) 2 (lTeV/m ) 3 sec, where £ < 1 is a 
numerical factor || . For relevant values of r\ and this 



decay is essentially 'instantaneous' on cosmological time 
scales at late epochs of interest to us. 

By far the largest number of particles eventually pro- 
duced by an X would come through the hadronic jet frag- 
mentation of quarks and gluons resulting from its decay 
(see, e.g., Ref. |l8| for arguments concerning the domi- 
nance of the hadronic decay channel) [[l9| . The fragmen- 
tation of the quarks/gluons into jets of hadrons and the 
photon spectrum resulting from decay of neutral pions 
in these jets are well described by the string fragmenta- 
tion scheme implemented in the JETSET program [ p0| . 
We assume typical hadronic 3-body decays jl8j of a X 
particle into all kinematically available quark pairs and 
one uncolored (assumed massless) spectator, and obtain 
the injected photon spectrum from the decay of a sin- 
gle higgs X particle by using the parametrization of the 
photon spectrum derived from JETSET, as described in 
Ref. (HI . Folding this spectrum with the X particle pro- 
duction rate then gives us the full injection spectrum. 

In our calculation of the predicted total 7-ray flux to- 
day, we have included the effects of electromagnetic cas- 
cading and 7-7 scattering and also included the 
effect of absorption due to pair production on photons of 
infrared, optical and ultraviolet backgrounds (2^]. The 
details of these calculations will be given elsewhere. In 
our numerical calculations we have assumed a spatially 
flat universe with fig — 1 and present Hubble constant 
H = TSkmsec" 1 Mpc -1 . 
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FIG. 1. Gamma ray flux due to decay of the higgs from 
TD processes with p = 1, higgs mass = 500 GeV (dashed 
curve), 1 TeV (solid curve), and 2 TeV (dotted curve), and 
Qo,<p = 4.5 x I0~ 23 eVcm -3 sec -1 . The extragalactic dif- 
fuse 7-ray background data from EGRET (Ref. [IT]) are also 
shown for comparison. 
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Fig.l shows the higgs contribution to the EDGRB for 
TD processes with p = 1 and higgs mass = 500 GeV, 
1 TeV, and 2 TeV. The normalization of the curves in 
Fig. 1 corresponds to Qo,0 ~ 4.5 x 10~ 23 eV cm~ 3 sec -1 , 
which, as evident from Fig. 1, is an upper limit on Qo,0 
(and hence on Q ) for p = 1 TD processes imposed by 
the EGRET data. 

A recent analysis |l7| of the EDGRB indicates that 
the spectrum continues at least up to ~ 100 GeV. The 
EDGRB up to ~ 10 GeV can be interpreted as arising 
from a superposition of unresolved blazars p3[ . However, 
if blazar 7-rays are produced by Compton upscattering 
of lower energy blazar photons, then only X-ray selected 
BL Lac objects, a small fraction of the blazar population, 
may qualify as possible contributors to EDGRB above 10 
GeV, because only they may have the requisite luminosi- 
ties. The majority of radio selected BL Lac objects and 
flat spectrum radio quasars are likely to have luminosities 
falling off above ~ 10 GeV j24|. Also, extragalactic ab- 
sorption effects are likely to steepen the high energy spec- 
tra of high-redshift quasars above ~ 10 GeV |2^j . Thus, 
there may be cause to consider another component of cos- 
mic diffuse 7-ray emission. From the shapes of the curves 
in Fig.l, we see that decays of higgs of mass ~ 500 - 1000 
GeV from p=l TD processes in flat SUSY theories may 
play an important role in producing the EDGRB in the 
10 - 100 GeV energy range. The rate of energy injection 
in the form of TeV higgs needed to explain the EDGRB 
above a few GeV is Q EDGRB ~ 4.5 x 10~ 23 eV cm -3 sec" 1 
for m ~ 1 TeV. 

The upper limit on Qo,4> from EDGRB also implies 
(through energy equipartition arguments) an upper limit 
on Q as well as on Qo,v, i.e., Qo/2 ~ Qo,4> ~ Qo.y < 
4.5 x 10 -23 eVcm -3 sec -1 . There is, however, an inde- 
pendent upper limit imposed on Qo,v (and hence on 
Qo) by the HECR data: Q ,y < Qo ECR - 3.3 x 
10- 22 (? ? /10 16 GeV) a5 eVcm- 3 sec- 1 ||, where Q$ EGR 
is the rate of energy injection needed to explain the 
HECR. For 77 > 1.9 X 10 14 GeV, we overproduce EDGRB 
if we wish to explain the HECR, and is, therefore, unfa- 
vored. (Of course, TDs with 77 > 1.9 x 10 14 GeV can give 
significant contribution to the EDGRB while not con- 
tributing significantly to the HECR.) For 3 x 10 11 GeV < 
77 < 1.9 x 10 14 GeV we can produce the HECR, but 
in this case we significantly underproduce EDGRB if 
7/ <C 1.9 x 10 14 GeV. The two independent upper lim- 
its can be saturated, i.e., we can explain both EDGRB 
above a few GeV and HECR, if 77 « 1.9 x 10 14 GeV. 

The value of Qo for a general TD process is not known 
a priori — it depends on at least two (not necessarily 
mutually independent) unknown parameters, namely, the 
fraction of the total energy density of the relevant defects 
going into X particles, and the symmetry breaking scale 
77 at which the relevant TDs are formed. Therefore, the 
above arguments do not by themselves rule out the exis- 



tence of GUT scale (i.e., 77 ~ 10 16 GeV) TDs in general 
- they only tell us that GUT scale TDs are unlikely 
to be responsible for HECR, because that would conflict 
with EDGRB. 

The situation is, however, very different in the case 
of cosmic strings with the recent results pp| ] of numer- 
ical simulations of cosmic string evolution in the Uni- 
verse. These studies show that the energy density, p 8 (t), 
in "long" (i.e., horizon crossing) strings at any time t is 
maintained in the scaling solution [Q, p s (t) = /i/(x 2 t 2 ), 
by energy loss from long strings occurring predominantly 
on the scale of the string width, i.e. through formation 
of string width-size small loops which quickly decay into 
X particles or through direct emission of the X particles 
that 'constitute' the strings, and not through formation 
of (sub)horizon size loops and their subsequent decay by 
emission of gravitational radiation as thought earlier . 
(Here fj, is the energy per unit length of the string, and x 
is in the range 0.27 - 0.34 Jl(]].) This result, while subject 
to confirmation by independent simulations, obviously 
has important implications for HECR and EDGRB. In- 
deed, in this case, there is effectively only one free param- 
eter (namely, 77 or equivalently, 77) in the problem, which 
also fixes Qo- I n fact, in this case, the observed data 
on ultrahigh energy cosmic rays already rule out [p|JTc|] 
GUT scale cosmic strings for the standard non-flat po- 
tential case. We shall see that this is also true for the flat 
potential case, but here an additional constraint (due to 
TeV mass scale higgs) comes from the EDGRB data. 

From the results of Ref. |Q, the rate of energy loss 
of strings per unit volume through X particle emission 
is dpx/dt ~ (2/3)p/(x 2 t 3 ), which gives Q ~ 7.4/j/t; 3 , 
with x ~ 0.3. Requiring Q > ~ 0.5Q < Q EDGRB ^ 



4.5 x 10" 23 eVcm 



we get n < 4.5 x 10 



(10 16 GeV) 2 . Taking, for flat potentials, fi ~ O.I77 2 [|, 
we get 77 < 2.1 x 10 14 GeV. Thus, in this case, GUT scale 
cosmic strings necessarily overproduce EDGRB, and are, 
therefore, ruled out. A similar constraint follows from 
HECR: Here one requires Q ,v « 0.5Q < Qo ECR - 
3.3 x 10" 22 (77/10 16 GeV)°- 5 eV cm- 3 sec~ 1 , which gives 
77 < 2.2 x 10 14 GeV. 

Note that for a non-flat potential, where p ~ 7T77 2 and 
where one expects the X particles to be predominantly of 
heavy mass scale ~ 77, so that Qo.v ~ Qoi the constraint 
on 77 from HECR is 77 < 1.4 x 10 13 GeV. 

It is thus clear that cosmic strings with 77 much greater 
than 10 14 GeV, and in particular, GUT scale cosmic 
strings with rj ~ 10 16 GeV, are ruled out both for flat as 
well as non-flat potentials i/X particle production is their 
dominant energy loss mechanism. At the same time, in 
this case, cosmic strings with ?)~2x 10 14 GeV in SUSY 
models with flat potentials can potentially account for 
the high energy ends of both EDGRB and HECR. In 
this respect, absence of free parameters other than the 
symmetry breaking scale 77 seems to make cosmic strings 
a "natural" candidate source of HECR (and possibly of 
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EDGRB above a few GeV). 

Cosmic string formation at r\ ~ 10 GeV rather than 
at the GUT scale of ~ 10 16 GeV is not hard to imple- 
ment. For example, in a SUSY theory, the breaking 
SO(10) -> SU(3) x SU(2) x U(1) Y X U(l) can take place 
at the GUT unification scale Mqut ~ 10 16 GeV without 
any cosmic string formation, but the second U(l) can be 
subsequently broken with a flat potential with a VEV 
77 ~ 10 14 GeV to yield cosmic strings that are relevant 
for EDGRB and HECR. 

Cosmic strings with 77 ~ 10 14 GeV would be too light to 
be relevant for structure formation in the Universe. Their 
signatures on the CMBR sky would also be too weak to 
detect. However, signatures of these cosmic strings may 
be searched for with next generation 7-ray instruments 
such as GLAST [^6|, which will be able to resolve fur- 
ther the discrete source component and thereby reveal 
the possible existence of a truly cosmic component of the 
EDGRB as in Fig. 1, and in proposed HECR observato- 
ries such as Auger §7| and OWL §8|. 
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